New Branch-and-Cut Algorithm for Bilevel Linear Programming

Linear mixed 0–1 integer programming problems may be reformulated as equivalent continuous bilevel linear programming (BLP) problems. We exploit these equivalences to transpose the concept of mixed 0–1 Gomory cuts to BLP. The first phase of our new algorithm generates Gomory-like cuts. The second phase consists of a branch-and-bound procedure to ensure finite termination with a global optimal solution. Different features of the algorithm, in particular, the cut selection and branching criteria are studied in details. We propose also a set of algorithmic tests and procedures to improve the method. Finally, we illustrate the performance through numerical experiments. Our algorithm outperforms pure branch-and-bound when tested on a series of randomly generated problems. Work of the authors was partially supported by FCAR, MITACS and NSERC grants.     

Shock wave diffraction over a 90 degree sharp corner — Posters presented at 18th ISSW

A special poster session was run during the 18th International Symposium on Shock Waves, held on July 21 – 26, 1991, in Sendai, Japan. The purpose of this session was to compare various CFD schemes which are useful for simulating shock wave phenomena. A 2-D planar shock wave diffraction over a 90 degree sharp corner was selected as a bench mark problem. The problem specification and the output format are described below. Experimental outputs for the close flow conditions to the computational ones were also called for. Twenty six poster presentations were made at this session and the majority of the posters reflected a standard of excellence that warrants publication in the Shock Waves journal. Therefore these results are reproduced here. However, since there is a limit on the available page numbers, it was not possible to publish all the posters that were presented at the special poster session. We express our sincere thanks to all the participants in this session.     

Small silver particles in glass surface layers produced by sodium-silver ion exchange — their concentration and size depth profile

Small spherical silver particles in a surface layer of commercial flat glass were produced by means of sodium-silver ion exchange. In each volume element of the layer there is a Gaussian distribution of the particle diameters. The mean diameter increases with penetration depth. Within one individual sample it can vary from 4.5 nm immediately at the glass surface up to more than 50 nm at the end of the layer. Due to a special preparation technique the results were gained by microspectrophotometric measurements as well as by investigations carried out with the transmission electron microscope and the electron-probe microanalyzer on one and the same sample always as function of the penetration depth.     

Characterization and Thermal Behaviour of Garnets from Almandine–Pyrope Series at 1200°C

Hyperfine Interactions - The natural garnets from almandine (Fe3Al2Si3O12)–pyrope (Mg3Al2Si3O12) series with the iron to magnesium atomic ratio ranging from 0.2 to 1 were characterised and...     

Dislocation unpinning model of acoustic emission from alkali halide crystals

The present paper reports the dislocation unpinning model of acoustic emission (AE) from alkali halide crystals. Equations are derived for the strain dependence of the transient AE pulse rate, peak value of the AE pulse rate and the total number of AE pulse emitted. It is found that the AE pulse rate should be maximum for a particular strain of the crystals. The peak value of the AE pulse rate should depend on the volume and strain rate of the crystals, and also on the pinning time of dislocations. Since the pinning time of dislocations decreases with increasing strain rate, the AE pulse rate should be weakly dependent on the strain rate of the crystals. The total number of AE should increase linearly with deformation and then it should attain a saturation value for the large deformation. By measuring the strain dependence of the AE pulse rate at a fixed strain rate, the time constantτ _s for surface annihilation of dislocations and the pinning timeτ _p of the dislocations can be determined. A good agreement is found between the theoretical and experimental results related to the AE from alkali halide crystals.     

Synthesis and adsorption properties of gold nanoparticles within pores of surface‐functional porous polymer microspheres

Mesoporous polymer microspheres with gold (Au) nanoparticles inside their pores were prepared considering their surface functionality and porosity. The Au/polymer composite microspheres prepared were characterized by transmission electron microscope (TEM), X‐ray diffraction (XRD), and Brunauer–Emmett–Teller (BET) techniques. The results showed that the adsorption of Au nanoparticles could be increased by imparting the pore structure and surface‐functional groups into the supporting polymer microspheres (in this study, poly (ethylene glycol dimethacrylate‐co‐acrylonitrile) and poly (EGDMA‐co‐AN) system). Above all, from this study, it was established that the porosity of the polymer microspheres is the most important factor that determines the distribution and adsorption amount of face‐centered cubic (fcc) Au nanoparticles in the final products. Our study showed that the continuous adsorption of Au nanoparticles with the aid of the large surface area and surface interaction sites formed more favorably the Au/polymer composite microspheres. The BET measurements of Au/poly(EGDMA‐co‐AN) composite microspheres reveals that the adsorption of Au nanoparticles into the pores kept the pore structure intact and made it more porous. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5627–5635, 2004     

Electrochemical behavior of Li intercalation processes into a Li[NixLi(1/3−2x/3)Mn(2/3−x/3)]O2 cathode

Li[Ni_xLi_(1/3−2x/3)Mn_(2/3−x/3)]O₂ (X=0.17, 0.25, 0.33, 0.5) compounds are prepared by a simple combustion method. The Rietvelt analysis shows that these compounds could be classified as having the α-NaFeO₂ structure. The initial charge-discharge and irreversible capacity increases with the decrease of x in Li[Ni_xLi_(1/3−2x/3)Mn_(2/3−x/3)]O₂. Indeed, Li[Ni_0.50Mn_0.50]O₂ compound shows relatively low initial discharge capacity of 200 mAh/g and large capacity loss during cycling, with Li[Ni_0.17Li_0.22Mn_0.61]O₂ and Li[Ni_0.25Li_0.17Mn₀.₅₈]O₂ compounds exhibit high initial discharge capacity over 245 mAh/g and stable cycle performance in the voltage range of 4.8 -2.0 V. On the other hand, XANES analysis shows that the oxidation state of Ni ion reversibly changes between Ni²⁺ and about Ni³⁺, while the oxidation state of Mn ion sustains Mn⁴⁺ during charge-discharge process. This result does not agree with the previously reported ‘electrochemistry model’ of Li[Ni_xLi_(1/3−2x/3)Mn_(2/3−x/3)]O₂, in which Ni ion changes between Ni²⁺ and NI⁴⁺. Based on these results, we modified oxidation-state change of Mn and Ni ion during charge-discharge process.     

A column generation approach for the timetabling problem of Greek high schools

A column generation (CG) approach for the solution of timetabling problems is presented. This methodology could be used for various instances of the timetabling problem, although in this paper the solution of the high-school situation in Greece is presented. The results obtained show clearly that the CG approach that has been extremely successful in recent years in the solution of airline crew scheduling problems could also be very efficient and robust for the solution of timetabling problems. Several large timetabling problems corresponding to real problems have been successfully solved, with the solutions obtained feasible and of very high quality in accordance with the problem definition. In addition, none of the solutions contained any idle hour for any of the teachers, which was one of the main goals of this optimization effort.     

20-relative neighborhood graphs are hamiltonian

For any two points p and q in the Euclidean plane, define LUN_pq = { v | v ∈ R², d_pv < d_pq and d_qv < d_pq}, where d_uv is the Euclidean distance between two points u and v . Given a set of points V in the plane, let LUN_pq(V) = V ∩ LUN_pq. Toussaint defined the relative neighborhood graph of V, denoted by RNG(V) or simply RNG, to be the undirected graph with vertices V such that for each pair p,q ∈ V, (p,q) is an edge of RNG(V) if and only if LUN_pq (V) = ?. The relative neighborhood graph has several applications in pattern recognition that have been studied by Toussaint. We shall generalize the idea of RNG to define the k-relative neighborhood graph of V, denoted by kRNG(V) or simply kRNG, to be the undirected graph with vertices V such that for each pair p,q ∈ V, (p,q) is an edge of kRNG(V) if and only if | LUN_pq(V) | < k, for some fixed positive number k. It can be shown that the number of edges of a kRNG is less than O(kn). Also, a kRNG can be constructed in O(kn₂) time. Let E_c = {e_pq| p ∈ V and q ∈ V}. Then G_c = (V,E_c) is a complete graph. For any subset F of E_c, define the maximum distance of F as max_epq∈Fd_pq. A Euclidean bottleneck Hamiltonian cycle is a Hamiltonian cycle in graph G_c whose maximum distance is the minimum among all Hamiltonian cycles in graph G_c. We shall prove that there exists a Euclidean bottleneck Hamiltonian cycle which is a subgraph of 20RNG(V). Hence, 20RNGs are Hamiltonian.